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OO Abstract 

1^ 

Q^ In the NMSSM, a light CP~odd Higgs arises due to spontaneous breaking of approxi- 

CO mate symmetries such as Peccei-Quinn or R-symmetry and is motivated by string theory. 

l/^ The case when it is heavier than two muons is well studied and constrained. We analyze 

^^ various meson decay, g — 2, beam dump and reactor bounds on the CP-odd Higgs with 

^^ mass below the muon threshold, in particular, addressing the question how light a CP- 

. . odd Higgs can be. We find that it has to be heavier than 210 MeV or have couplings to 

^ fcrmions 4 orders of magnitude below those of the Standard Model Higgs. Our analysis 

k>< applies more generally to couplings of a light pseudoscalar to matter. 
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1 Introduction 

An attractive extension of the minimal supersymmetric Standard Model (MSSM) is the next- 
to-minimal supersymmetric Standard Model (NMSSM). Historically, it has been motivated 
by the //-problem of the MSSM which can be solved by replacing the /x-parameter with an 
SM singlet S [1]. More recently, it has been suggested that the fine-tuning problem of the 
MSSM can be alleviated or removed in the NMSSM when a light CP-odd Higgs is present 
in the spectrum [21 [3] . Although this scenario is now tightly constrained by the new ALEPH 
analysis [Ij as well as BaBar data on T(3S') decays [5l ]6], some parameter space remains 
available [7] and can further be probed by ry^-decays [^. 

The NMSSM-like models can also be obtained from the heterotic string |9j, based on the 
MSSM constructions of Refs. [TOl [Til [12l [13] . An interesting feature of these models is that the 
SM singlet S is not a singlet under the full gauge symmetry of the heterotic string [E% x E%). 
Therefore all self-interactions 

are forbidden by gauge invariance. They are only allowed after spontaneous breaking of gauge 
symmetry and the corresponding couplings are suppressed by VEVs of the symmetry-breaking 
fields. On the other hand, interactions like SH1H2 can be allowed without symmetry breaking. 



Due to this hierarchy in the couphngs, one obtains specific versions of the NMSSM such as 
the Peccei-Quinn (PQ) or decouphng versions [9j. In the former case, a pseudo-Goldstone 
boson appears in the spectrum at the electroweak scale. Its mass is generated by small PQ 
violating effects and can be much below the GeV scale. 

Motivated by these considerations and also by a possible connection of a light pseudoscalar 
to dark matter [lU HSl [161 [T7] , iii this work we study constraints on a very light CP-odd Higgs 
and address the question how light a CP-odd Higgs can be. We analyze bounds from various 
meson decays, muon g — 2, beam dump and reactor experiments for the pseudoscalar mass 
below the muon threshold. Our results apply beyond the NMSSM to couplings of a light 
pseudoscalar to matter. 

2 NMSSM and a light CP-odd Higgs 

The NMSSM is the MSSM extended by a singlet superfield S. In what follows, we will focus 
on a particular version of the NMSSM which has no direct //-term, the so called Zs-symmetric 
NMSSM. The relevant superpotential of the Za-symmetric NMSSM is 

W = XSH1H2 + -kS^, (2) 

while the soft terms are given by 

1/soft = ml\Hi\^ + ml\H2\^ + m||5|2 + [\AxSH1H2 + '^kA^S^ + h.c.) . (3) 

A light pseudoscalar A^ appears naturally in two limiting cases: when the Higgs potential 
possesses either approximate Peccei-Quinn (PQ) or approximate R-symmetry [I8l [I9l EOj 1211 
El [211 [23]. 

2.1 Peccei-Quinn limit 

In the limit k — t- 0, the Lagrangian is invariant under the transformation 

i7i,2 ^ e*"Fi,2 , 5^e-2*"5. (4) 

At the electroweak scale this symmetry gets broken, resulting in the "axion" ^24j 

A"^ = ^(y sin2/3<sSM-2s5,) , 



N = yr;2sin2 2^ + 4s2 , (5) 

where ^^SSM — cos/3 Hn + sin/3 H21 is the MSSM pseudoscalar, s = {S) and the subscript 
I refers to the imaginary part of the Higgs neutral component. As usual, tan/3 = v\jv2 and 
V = \Jv\ + v\ = 174 GeV. The mass of the pseudoscalar is most easily expressed in the large 
tan^S regime [TO] : 

m^o — -'ii^Af.s . (6) 



Since renormalization of k is proportional to k itself, this coupling can be very small and 
A^ very light. In the string NMSSM example of Ref. [9J, k < 0{{4>/Mpi)^) with (/> being an 
average VEV of certain SM singlets. For (f) an order of magnitude below the Planck scale, k 
can be as small as 10~^ leading to a 100 MeV pseudoscalar. In other models it can be even 
lighter. Since the PQ symmetry is anomalous as in the DFSZ construction [25], the lower 
limit on m^o is set by the anomaly contribution of order 100 keV [25J (for s ~ v). 

2.2 R— symmetry limit 

In the limit A^^, Ax — )• 0, the Higgs sector of the NMSSM is R-invariant. Under R-symmetry 
the superfields transform as 

i/i,2 ^ e^°«//i,2 , S^e^°«5, (7) 

such that the superpotential transforms with charge 2. Spontaneous breaking of this symmetry 
results in an "R-axion" . Its composition is given by 



^° = ^(« sm2f3 AlissM + sSi') , 



N = y t>2 sin2 2/3 + s2 , (8) 

with ^SlsSM as in Eq. [s) 

Unlike the Peccei-Quinn symmetry, R-symmetry is not a (classical) symmetry of the full 
Lagrangian. Even if Af^,Ax — )• 0, the gaugino mass terms break it explicitly. Non-zero A- 
terms are induced by renormalization, so their minimal value is a loop factor times the gaugino 
mass. The axion mass is again approximated by ([6|. 

In both PQ- and R-symmetric cases, the light pseudoscalar is mostly a singlet in the 
limit s 3> u sin 2/3. Its couplings to gauge bosons and SM matter are suppressed in this limit. 
However, s cannot be too large, otherwise a large effective //-term is induced. An exception is 
the case A <C 1, which corresponds to the "decoupling limit", i.e. no communication between 
the singlet and the rest of the NMSSM. 

3 Constraints on a light CP— odd Higgs 

Following the notation of Ref. |7], the coupling of the CP-odd Higgs A^ to fermions is given 
by 

AC = -i- CAff ( rud d-i^d + - — 2^m„ u-y^u + mi l-f^l ] A^. (9) 

Zm\Y y tan p / 

In the NMSSM, the coupling CAff can be expressed in terms of the singlet-doublet mixing 
angle 6a and tan/3: Ca// = cos ^a tan /3 [7], with cos^^i = v sin 2/3 /A^ and N given in 
Eqs. |5|8[ In what follows, we treat it as a free parameter and derive various particle physics 
constraints on it. In the NMSSM, very large (> 10^) and very small (< 10~^) values of CAff 
lead to violation of perturbativity and/or finetuning, so it usually suffices to focus on the 



moderate CAff window. However, our analysis applies more generally to the coupling of any 
pseudoscalar to matter as long as the C^// parameter is universal for all fermions. 

For rrij^o > 2m ^, the resulting constraints, in particular from meson decays, are well 
studied. In the range 2m^ < rrij^o < 3m,r) ^^ decays predominantly into 2 muons, which 
is constrained by K^ — )• tt~^A^ and B ^ K A^. The resulting bound is of order CAff < 
0(10"^) [26j . Above the 3-pion threshold, the branching ratio for the decays into muons 
reduces and the bound weakens somewhat. For even larger A^ mass, the T — )• 'jA^ decay 
imposes CAff < 0.5 [7] at tan/3 ~ 1 with the bound getting weaker, 0(1), close to mr- Above 
12 GeV, the DELPHI data on e+e" -^ blA^ -^ bbbb set a hmit Ca// < 0(10) [7J. Further 
constraints, usually relevant at large tan/3, are summarized in Ref. [24J. 

The niAO < 2m^ territory is less well explored. Some constraints have been studied in 
Refs. [SgEO] in the framework of the NMSSM, and in Refs. [ZZIIMI for 2 Higgs Doublet Models 
(see also |29]). In what follows, we delineate the {mAO,CAff} parameter space taking into 
account meson decay, muon g — 2, nuclear reactor and beam dump constraints. In particular, 
we address the question how light a CP-odd Higgs boson can be. Since we work in terms of 
the coupling CAff, most of our results are largely independent of tan/3. 

For ruAO below the muon threshold, A^ can only decay into electron-positron pairs and 
photons. Its total decay width is 



r(A° ^ e+e-) + r{A° -^ 77) (10) 



with 




m'^ff) = ^^}-^Aoci,,^i-4-f, (11) 

r(A0^77) = "^^^f m\o \Y,rCAnQlhFiki) ' , (12) 

i 

where r = l{Nc) for leptons (quarks), ki = mf/m?.o, Qi is the charge of the fermion in the 
loop; Cau = CAff for the down-type fermions and Cau = CAff/tsLU^P for the up-type 
fermions. (Here we neglect the chargino contribution as the coupling is dominated by the SM 
fermions.) The loop function F{ki) is given by |30j 



-2(arcsin-^)2 for ki > I 



4' 



nk.) = { in„.i+yi^; ,„_.2 ,„„..! (13) 



5HSf^) + -P f-^ 



< 



and has the limits 



vT^m^^""^ iwxa, ^4, 



for ki < 1, 

,2 



kiF{ki)={ -^ forh = l (14) 



8 ^^'' '"« — 4' 
'2 



for ki > 1. 



For example, with m^o = 0.5 MeV and Ca// = 1 the total width is about 4 x 10~^^ eV. 
This corresponds to the decay length (for a boost factor 7 ~ 1) re ~ 60 km. Above the 
electron threshold, taking tuao = 50 MeV and Ca// = 1, the total width is 10~^ eV and the 
corresponding decay length is re ~ 2 cm. 
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Figure 1: tan/3 dependence of BR(^'^ — )• e+e ). The lowest curve corresponds to tan/3 = 1, 
the higher one to tan/3 = 3, and the highest curve to tan/3 = 10. 



The tan/3 dependence of BR(j4'^ — )• e~^e~) is shown in Fig. ^ By increasing tan/3 one 
reduces the up-type quark contributions to r(^'^ — )• 77), thereby increasing BR(A'^ — )• e+e"). 
This dependence saturates for tan/3 > 3. The decay mode A^ — )• e~^e~ dominates for rrij^o 
below 80 MeV and switches off abruptly just above the electron threshold. 

Below we study various constraints on {rrij^o, C^ff}- An important class of constraints is 
due to the decays 

X — )■ y + invisible 

with X, Y being some mesons. Experimental limits on their branching ratios exclude parts 
of parameter space where A^ is sufficiently long-lived to escape the detector, that is t^c > d 
with r being the lifetime oi A^, 7 a boost factor and d the size of the detector (~ 10 m). In 
other words, 



tot 



< 



rrij^o d 



For a two-body decay X ^ Y + A^, 



Eao 



m\ — m\ + m^o 
2mx 



and in relevant cases mj^ » my, ttt-J^o- Then for m^o < 2me, one has 



"^A<'V<^A// < 25 MeV X 



mx 
GeV 



(15) 



(16) 



(17) 



for d ~ 10 m and tan/3 ~ 1. For rrij^o > 2me, one can estimate the resulting bound by taking 
Ftot ~ Tg+e-, in which case 



rriAoCAff < 8 MeV x 



mx 
GeV 



(18) 



3.1 Rare B— decays B -^ K + invisible 

Limits on production of A^ in rare i?-meson decays [311 [32l [33l [M! result from the following 
bounds on the branching ratios measured by CLEO f35] and BaBar j36j : 

^CLEO^^o ^ kO ^ invisible) < 5.3 x 10"^ , (19) 

^BaBar^^- ^ K-ui?) < 7.0 X 10"^ , 

where in what follows we use the more constraining CLEO result. Since in both experiments 
A^ appears as missing energy, the bounds apply only if A^ decays outside the detector. For 
mx = m^o = 5.28 GeV, this implies 



(■e ) 



nij^o \JCAjf ;$ 37 MeV for tti^o < 2me 
m^o Ca// < 18 MeV for m^o > 2m^ . (20) 

In the NMSSM, the decay rate for 5° -^ K^A^ is given by 



r(i?o -. K^A^) = "^;"^' \CA' ^ /o-"(mio) p^i^^y, (21) 



210^5 i^-i ^1^ 



2 2 



mf, 



where the form factor /q (0) ~ 0.3 — 0.4 [37j and \pk\ — m.^o/2 is the three momentum 
of the kaon. The quantity Ca has been calculated in Ref. |26j in the large tan/3 regime, 
Ca ~ ^Aff tan /3 mhijit for order one stop mixing and EW scale sparticles. Since the full 
NMSSM calculation at low tan /3 is not available, we estimate the order of magnitude of the 
resulting bound by a rescaling of this result. Using the total B^ width F^o = 4.3 x 10~^^ GeV, 
the CLEO bound implies Caj/ < 0.02/ tan /3. Taking conservatively tan/3 ~ 0(1), we geir] 

CAff < 10-2 . (22) 

This constraint is already strong at small tan/3 and gets even stronger at large tan/3. 

The resulting exclusion region is shown in Fig. ^ (marked "i?*^ — )• K^ + inv."). We note 
that, in contrast to the lower boundary, the right boundary of this region is calculated quite 



reliably from Eq. 15 and is essentially independent of tan/3. In the plot, we use the full A^- 
width Ftot without resorting to the approximation Ftot ~ Te+e-- The kink at rra^o « 2me is 
due to the rapid fall of T{A^ — )• e~^e~) as niAO approaches the threshold from above. Finally, 
the dependence on the detector size is only square-root. 

3.2 Rare K— decays A' — )■ vr + invisible 

A light (invisible) A^ can also be produced in i^-decays. The relevant branching ratio has 
been measured by E787 j38l [39] and E949 [4OJ0 

^E787^^+ ^ 7r+ + invisible) < 4.5 x 10-^\ (23) 

i3^9^^(K+^7r++ invisible) < 10"^°. 



^Essentially, this corresponds to the SM contribution with an additional coupling ([9|. The bn — SL transition 
is mediated by the W — t loop with A^ coupled to the top quark. 



^These bounds become significantly weaker at the pion pole, m^o = m^ 



We will use the tighter E787 bound. Eqs. 17 and 18 for mx = rnK+ — 494 MeV yield the 
"invisibility" conditions 



"-e ) 



ruj^o \/CAff ;$ 21 MeV for tti^o < 2me 
mjifi CAff < 5 MeV for m^o > 2me ■ (24) 



The decay rate is given by 



g'f \vtsv;j' ,^, ,2 |p 



r(K+ ^ vr+A") = ^^' ' ;i *'^' \C'a\'' -^ /o^^Ko) R^^— ^ , (25) 



2i07r5 ■ --■ ^^_^ 



fA'+, 



2 /m^_|_ — ?n. 



m.. 



where the form factor f^ (0) ~ 1 [H], IpttI — m^+/2 and C^ ~ Cyi//tan/3 rngmt. With 
r^+ = 5.32 X 10~^^ GeV, the experimental bound requires Caj/ < 2 x 10~^/tan/3, which for 
tan/3 ~ 0(1) gives 

C^// < 10-^ . (26) 

The corresponding excluded region is marked "X^ — t- vr^ + inv." in Fig. [2] As in the case 
of i?-decays, this bound only gets stronger with increasing tan /? and its precise value is not 
important for us. 

3.3 Rare decays B —^ K e^e^, fC — )■ tt e+e^ 

If A^ decays inside the detector, it contributes to the processes B ^ K e~^e~ , K ^ it e^e~ . 
In experimental measurements, one sets a cut on the invariant mass rrig+g- > 140 MeV in 
order to suppress backgrounds including conversion photons and tt'^ — >• e^e~7. Therefore, the 
resulting bounds apply for tti^o > 140 MeV. 
BELLE has reported 



^BELLE^^ ^ ^ ^+^-) ^ (4.8+J:° ± 0.3 ± 0.1) X 10-^ , (27) 

where i includes muons and electrons with m^+f,- > 140 MeV. Assuming lepton universality, 
we win use a conservative bound BR{B — > K e+e') < 2.4 x 10"'''. The NMSSM result can 



be read off from Eq. 21 using 

BR(S ^ K e+e") ~ BR(5 -^ K A^) ^ BR(A° -^ e+e") . (28) 

In the relevant mass range, BR(yl'^ — )• e^e") is 20 - 40 % for tan/3 ~ 1. After imposing the 
condition that A^ decay inside the detector, we get 

CAff < 8 X 10"^ , (29) 

which excludes the strip 140 MeV < my^o < 2m^ shown in Fig. ^ (marked "S — t- Ke^e~^^). 
A similar result is obtained from X-decays [43j: 
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Figure 2: Constraints from meson decays and muon g — 1. The colored regions are excluded. 
These bounds include the effect of varying BR(yl'^ — > e+e") with 777.^0. 



which also employs the same kinematic cut rn-g+g- > 140 MeV. Using Eq. 25 and an analog 



of (28), we get 



CAjf < 2 X 10-2 , (31) 

where the "visibility" condition for A^ has been imposed. This bound is shown in Fig. [2] 
marked "i^ — >• 7re''~e~". 

It is noteworthy that the window 140 MeV < m^o < 2m^ is eliminated by 2 different 
processes. There are additional NMSSM contributions to B ^ K e^e~, K ^ tt e~^e~ apart 
from that of A^, so in principle there could be cancellations. Considering two independent 
processes makes this possibility less likely. 

The BaBar measurement of BR(i? — )• K i~^l~) [H] imposes a lower kinematic cut 
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?7ig+g- > 30 MeV thereby losing somewhat in efficiency because of the low energy back- 
grounds [26j . Their result 

^BaBar^^ ^ ^ ^+^-^ ^ ^Q_34 ^ Q Q7 ^ Q Q3^ ^ ^Q-6 (32) 

excludes the region m^o > 30 MeV with Ca/j > 10~^ — 1 (where A^ decays inside the 
detector) depending on the pseudoscalar mass, Fig. ^^ (marked "-B — )• Ke~^e~"). Although one 
may question the reliability of this result at low e~^e~ invariant masses, another experiment, 
to be discussed in the next subsection, excludes a similar region of parameter space. 

3.4 Rare K-decays K+ -^ 11+ + X 

A byproduct of the i^^2 experiment in Japan was a measurement of a 2-body decay K~^ — )• 
vr"^ + X, where X is any particle [35]. One searched for a peak in the 7r+ momentum for 10 
MeV < mx < 300 MeV. The resulting bound is 

B{K+ ^ TT-^ + X) < IQ-^ (33) 

at 90% CL for mx < 60 MeV. For larger mx up to 120 MeV this bound relaxes to 10~^. 

The excluded parameter space is shown in Fig. p^ (marked "if+ — )• vr+X"). The constraint 
amounts approximately to 

CAff < 4 X 10-2 (34) 

for m^o > 10 MeV. 

We note that at m^o — tti^o the constraint is weaker. However, this region is disfavored 
by vr^ — )• e^e~ (see below) and the vr"^ — vr*^ mass difference when vr'^ — A^ mixing is taken into 
account. 

A similar region of parameter space (up to m^o = 100 MeV) is excluded by the process 
7r+ — )• e^u A^ with subsequent decay ^4'^ — )• e+e" [46J (for applications to axion models, 
see [17]). 

3.5 Radiative Upsilon— decays 

The bounds on A^ production in radiative T decays come from CLEO [48j and BaBar ^491 
EDI: 



^cleO(y(;L5) ^ 7 + invisible) < 1.3 x W'^, (35) 

gBABAR('^(3^) ^ 7 + invisible) < 3 x 10^^ (36) 

They apply only if A^ decays outside the detector, which for m,x = n^T(3S) — WA GeV means 

iriAO \/CAff ^ 44 MeV for m^o < 2?7ie , 
m,A0 CAff ;$ 25 MeV for tti^o > 2me . (37) 

The branching ratio for T — )■ A^^ is given by [51^ [52l [26] 

B{T ^ ^07) _ Gpml 2 f, ^' 

: ^Aff i i - - 



K/-nr . — 4~"^ ~ /TT ^Aff { 1 4" ) -^QCD , (38) 



10 



where Fqcd ~ 0.5 is a QCD correction factor, B{T{1S) -^ n+fi') = 0.025 and B{T{3S) -^ 
fi+lj,-) = 0.022 [53]. The resulting bounds (marked in Fig. [2] "T(15) -^ 7 + inv." and 
"T(35)^7 + inv.") are 



CAff < 0.37 (CLEO) , 
CAff < 0.19 (BaBar) , 



(39) 



independent of tan /3. 



3.6 Pion decay tt*^ — t- e+e 

A hght A^ provides a pseudoscalar channel for pion annihilation into e+e^ (see, e.g. [.54j). 
This chirality-suppressed decay proceeds in the SM through a loop diagram with a 7r'^77* 
vertex and has a very small branching ratio. The recent KTeV result |55j 



B 



KTeV/ 
vr 



e e 



(7.48 ±0.29 ±0.25) x 10 



-8 



(40) 



is somewhat {3a) above the SM prediction \56\ [57 ] . To be conservative, in what follows we 
will require that the tree level contribution from A^ not exceed the central experimental value. 



We find 



< 7.5 X 10"^ 



r(7r^ 



^0 



e+e-1 



Gp 



At: Iml. 



mimlf^ 



C\ 



m 



AO 



±ir^om^o|2 ^ff ' 



(41) 



where we have neglected the up-quark contribution, {Olrudd'y^ d\TT^) ~ — i?7i^/7rjj T^o is given 
by Eq. 10, m-j^ = 135 MeV and fi^ = 93 MeV. We neglect the vr'^ — A^ mixing effects which 
are of order 5m? /m^ ~ f-j^jM^^ ~ 10"^ and relevant only very close to the pion mass. 



The total width of vr*^ is 



r(7r° -^ 77) 



a 



m- 



647r3 /2 



then the KTeV result requires 



Cah < 20 



(42) 



(43) 



away from the resonance region, where the constraint is stronger. A more precise bound 
including the m^o-dependence is shown in Fig. [2| marked "vr'' — )• e'^e~". This constraint 
is complementary to those of the X ^ Y + invisible decays, in that it excludes parameter 
space above CAff — 20 regardless of the ^4*^ mass. It is also a reliable tree-level constraint 
essentially independent of tan (3. 



■^Our bound on Ca// is not sensitive to this approximation as it scales as a square root of this matrix 
element. 
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3.7 Muon anomalous magnetic moment 

At loop level, A^ contributes to the muon g — 2 which is well measured. Currently, there is 
a 4(7 discrepancy between the SM prediction for the muon g — 2 and its measured value at 
BNL E821 [58]: 

Aa^, = oj^^P - a™ = (31.6 ± 7.9) x 10-^° , (44) 

which may be considered a hint for new physics. 

In the NMSSM, there are significant one- and two-loop contributions of the CP-odd Higgs 
A^ to a^. They are given, for example, in |59j : 



6a^{A^) = 6a]!^{A^) + 5a'^{A^), (45) 

aL//iO\ _ V2Gf _2 ir^ |2 f f^Ao' 



V2Gf 2 i^ |2 . fn^^ 
-^-^- ^1, \CAff\ Jl — 2 



K(^) = -''^<\CAffVh 



V2GFa 



5af{A^) = ^-^^ml\C,~^- 



1 „ / ml \ 1 J, / fn^ \ „ / w-? 



Q+ 2 «-^2( ^2~) + Q'^2( ^2~) +-^^\ 2 

3 tan^ p Vm^o/ 6 \m.Q/ \m 



where 



fi{z) = / dx^^ — ^ , 

Jq x^ + z{l-x) 

f2{z) = z f dx .^ \ in^ii^. (46) 







x{l — x) — z 



The one loop-contribution is negative which makes the discrepancy worse. For tti^o above 
roughly 1 GeV, the two-loop contribution may be dominant and resolve the discrepancy. 
However, this does not occur in the mass range we consider. 

Since there are NMSSM contributions to (7 — 2 of both signs, the contribution from the 
CP-odd Higgs A^ can be canceled. We then require that the latter not worsen the discrepancy 
beyond 5 a: 

5a^,{A^) < a|^"P - af^ ~ 40 x 10-^° {5a) . (47) 



The corresponding bound on C^// (at tan/3 r^ 1) is shown in Fig. l2l marked "a^". It can be 
approximated by 

CAff < 2 (48) 

for m^o < rn,^. The tan/3 dependence is very mild in the region of interest and stems only 
from the 2-loop contribution, which is subdominant. 

Once this bound is imposed, the electron g — 2 constraint is satisfied automatically. 

3.8 Other constraints 

Further (model-dependent) constraints are summarized in Refs. |53l[26]. These are weaker 
than the bounds we have considered and require assumptions about the NMSSM spectrum. 
For instance, there are contributions from all neutral Higgses to Bs — )• /U^/i^ and B — B 
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mixing which ahows one to eUminate parts of parameter space {Ca/j ~ 0(10)) depending on 
their masses and tan/3 [26j . 

Among other possible flavor physics constraints are J/^ decays. CLEO has recently 
reported B [J/"^ — )• 7 + invisible) < 4.3 x 10~^ [60], which is somewhat weaker than the 

analogous T(35) bomid. The ^4" couphng to up-type quarks falls very quickly with tan/3, so 
we do not use this result in our analysis. 

Further, the missing-energy process B — )• KA^A^ proceeding through the hA^A^ cou- 
pling [HI] sets a mild constraint on the SH1H2 coupling in the superpotential, A < 0.7. 

A light CP-odd Higgs could potentially be constrained by the LEP data. ^4*^ couples to 
the Z-boson at tree level through the A^Hj^Z^ vertex |24j, therefore the (invisible) Z-width 
does not constrain the mass of A^. Furthermore, the electroweak oblique corrections are 
suppressed by the mass of the heavier pseudoscalar (see, e.g. [IS]). For the same reason, the 
A^ production at LEP through e^e~ -^ h A^ is suppressed. The associated production with 
bottom quarks e~^e~ — )• bb A^ is also insignificant [15j. Finally, the constraints from Z — )• 7 ^4'^ 
are weak [62l [63] . 

Astrophysical bounds have been summarized in Ref. [20]. They are usually relevant for 
sub-MeV pseudoscalar masses which, in the range 10~^ < Cajj < 10^, are already excluded 
by meson decays with missing energy. However, the supernova SN1987A sets an additional 
constraint for a small coupling: Ca// > 10~^ when tti^o < 30 MeV [20] . 

To summarize this section, we see that a combination of various constraints requires the 
CP-odd Higgs to be heavier than 2m^ (unless Ca// < 10~^). To obtain this bound we did 
not rely on the specifics of the NMSSM. All we used was the coupling ^ at tan/3 ~ 0(1), 
which is much more general. This coupling is sufficient to induce the b — s and s — d transitions 
(with flavor change due to SM loops) which were used in the processes like B ^ K A^ and 
i^ — 7- vr A^. Similarly, T decays, vr^ — >• e^e~ and muon g — 2 are generated directly by (|9|. 

For completeness, in the next section we discuss the reactor and beam dump results which 
have been used in the past to constrain axion models. 

4 Further bounds from reactor and beam dump experiments 

4.1 Reactor bounds 

Searches for axion-like particles using nuclear power reactors set constraints on the parameter 
space of the CP-odd Higgs. Here we consider 2 representative experiments which employ 
Bugey and Kuo-Sheng nuclear reactors. 

Axion-like particles can be emitted in place of photons from excited nuclear levels which 
makes nuclear reactors a source of pseudoscalars with masses up to 10 MeV. In [6l], the 
detector was placed 18.5 m from the Bugey reactor core and one searched for the decays 
A^ — )• e~^e~ . No excess of e^e~ events has been observed which set a constraint on the axion 
decay constant /;^. The corresponding exclusion region can be read off from Fig. 5 of [6l] 
using the conversion 

_ 1 2mw ... 

^Aff - -7 — • (49) 

jx y 

The result is shown in Fig. [3] 

13 



10^ 

102 

10 = 



d 



10 



10- 



10" 



- Kuo-Sheng 

1 1 1 1 1 1 1 


\ 


Bugey 

1 1 1 1 1 


1 1 1 1 1 1 1 1 



-1 ^ 



10" 



10"^ 10-2 

m^o [GeV] 



10-1 



Figure 3: Constraints from nuclear power reactors. 



Another experiment at the Kuo-Sheng nuclear reactor searched for axions via Compton 
conversion on electrons [65]. A Germanium detector, placed 28 m away from the reactor, 
measured the ionization energy resulting from the axion-photon conversion in the detector. 
Again, no signal was found. The exclusion region can be read off from Fig. 7 of [65J and their 
Eq. 31, gaee qInn < 1-3 X 10~i°. Using gaee = CAffgme/{2mw) and g\j^j^ 
this translates into 



niAo Cajj < 2 X 10"'' MeV 



3xlO-*^m^o/eV, 
(50) 



for rriAo < 2me- The experiment is also sensitive to tti^o up to 2.23 MeV, which is the 
pn — )• d-f transition energy. Requiring that the axion not decay before it reaches the detector 
{niAoCAff < 0.3 MeV), one obtains the bulge at m^o > 2me in Fig. Isj 



4.2 Beam dump limits 

Axion-like particles can be emitted via bremsstrahlung or Primakoff production in beam 
dump experiments (see, e.g. [66]). The setup of these experiments is as follows. An intense 
beam of particles (electrons or protons) hits a thick target, which absorbs the beam and the 
interaction products apart from very weakly interacting particles such as axions. The decay 
products of the latter are collected by the detector, typically placed tens of meters behind the 
target. 
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Below we consider 4 representative beam dump experimentsj^] 

• SLAC E141 [68j: 2 x 10^^ electrons at energy 9 GeV struck a 12 cm tungsten target; 
detector 35 m behind the target 

• Fermilab E774 |69j : 0.52 x 10^'' electrons at 275 GeV dumped at a 30 cm target; detector 
length 7.25 m 

• CHARM [7D]: 2.4 x 10^^ protons at 400 GeV dumped at a thick copper target; detector 
480 m behind the target 

• Orsay |71j : 2 x 10^^ electrons at energy 1.6 GeV dumped in a 1 m target; detector 2 m 
behind the target 

The corresponding exclusion regions (Fig. H| can be read off from the plots presented in these 
papers either using the conversion factor for the axion decay constant ( 49 ) or calculating the 



axion decay time according to Eq. 10 




10-3 10-2 

m^o [GeV] 

Figure 4: Constraints from the beam dump experiments. 



The reactor and beam dump experiments by themselves eliminate most of the parameter 
space (Figs. 314). They are based on a different kind of physics compared to meson decays 



*We are not displaying the results of the SLAC E137 [SS] and Fermilab 605 [ST] experiments since the 
corresponding exclusion regions are largely covered by other experiments. 
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of Sec. [3] and in this sense are complementary. Note also that, unlike meson decays, electron 
beam dump experiments as well as the muon g — 2 probe directly the lepton-axion coupling, 
which could be the only axion coupling to matter in exotic ( "leptophilic" ) scenarios. 



5 Conclusion 

In this work, we have studied the question how light a CP-odd Higgs A^ of the NMSSM can 
be. We have analyzed constraints from meson decays, muon g — 2, beam dump and reactor 
experiments. We find that the parameter space m^o < 2m^ is excluded (by more than one 
experiment) unless the coupling of A^ to matter is 4 orders of magnitude smaller than that 
of the Standard Model Higgs, i.e. C^// < 10~^. Since such a small coupling can hardly be 
achieved in the NMSSM, we conclude that A^ has to be heavier than about 210 MeV. 

Our analysis applies more generally to couplings of a light pseudoscalar to matter. We 
have not used any specific features of the NMSSM (nor supersymmetry) . We have only relied 
on Eq. [9] and analyzed parameter space in terms of {m^o, Ca/j}- Flavor changing couplings 
of the Standard Model are sufficient to generate the processes like B ^ K A^ and K ^ tt A^, 
which lead to strong constraints. 

Since the CP-odd Higgs is heavier than 2m^, the decay channel to muons is open. One can 
therefore produce ^4'^ in gluon fusion at the LHC and look for fi^fj-^ pairs with low invariant 
mass [72] , or search for decays h — )• 2A^ — t- Afi |73| ^^ (similar to what has been done by DO 
at Tevatron [75]). This will provide an important test of models with a light pseudoscalar. 

Acknowledgments. A.R. would like to thank G. Weiglein for useful discussions. 
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